clear all;clc;close all;
c=3e8;
Rsat=25480e3;
Rl=6380e3;
tmax=43200;  %12 часов время оброта
w=2*pi/tmax;
Rmas=[nan;nan];
F0=1602e6;
Fdmas = 0;

for t=2:tmax*2
    Rsvx=Rsat*cos(w*t);
    Rsvy=Rsat*sin(w*t);
    R=[Rsvx;Rsvy];
    Rmas(:,t)=R;
    Rr=[0;Rl];
    Rv=[Rmas(:,t)-Rr];
    V(:,t)=(Rmas(:,t)-Rmas(:,t-1))/1;
    Vu=-dot(V(:,t),Rv)/norm(Rv); 
    Fd=(Vu*F0)/c;
    Fdmas(:,t)=Fd;
end
figure(1);
plot(Rmas(1,:),Rmas(2,:));
figure(2);
plot(1:tmax*2,Fdmas);
    

